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Abstract 

We investigate the time evolution of the decay (or ionization) probability 
of a D-dimensional model atom (D = 1,2,3) in the presence of a uniform 
(i.e., static and homogeneous) background field. The model atom consists 
in a non-relativistic point particle in the presence of a point-like attractive 
well. It is shown that the model exhibits infinitely many resonances leading 
to possible deviations from the naive exponential decay law of the non-decay 
(or survival) probability of the initial atomic quantum state. Almost stable 
states exist due to the presence of the attractive interaction, no matter how 
weak it is. Analytic estimates as well as numerical evaluation of the decay 
rates are explicitly given and discussed. 
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I. INTRODUCTION 



Exponential decay is a common feature of many physical processes; in particular, it is the 
universal hallmark of unstable systems such as radioactive nuclei. However, it is known that 
under very general conditions quantum mechanics predicts deviations from the exponential 
decay within short as well as long time intervals [1]. As pointed out by Khalfin [2], the latter 
situation occurs whenever the spectrum of the Hamiltonian H is bounded from below; in 
this case, the Paley- Wiener theorem [3] on Fourier transforms implies that the non-decay or 
survival amplitude A{t] [■?/']) := (■?/'(0)|'?/'(t)) necessarily satisfies 

This condition clearly rules out an exponential decay for i — > oo, as this would cause the 
integral above to diverge.^ Explicit calculations in a number of models [4] show that, in 
fact, there occurs a crossover from exponential to power law decay when t oo. 

One may wonder what happens if the Hamiltonian is not bounded from below. At 
first glance this might appear to be an academic question, since any realistic Hamiltonian 
should be bounded from below. However, such "unrealistic" Hamiltonians are often found 
in physics. Some examples are the decay of a metastablc vacuum through the formation 
of bubbles of the true vacuum [5], the droplet model for first order phase transitions in 
statistical physics [6], or the ionization of an atom by a static electric field [7]. In the latter 
case Herbst [8] provided a partial answer to that question. Let H — —A -\-V -\- Fx be 
the Hamiltonian describing a one-electron atom in a uniform electric field. If V{x, y, z) is 
holomorphic in x and V{x + ia,y,z) is bounded and decreases to zero as r := {x'^ + y'^ + 
^2^1/2 _^ ^ ^Qj, each, a e M, whereas -0 is an eigenvector of —A + V with negative eigenvalue, 
then (for F ^ 0) 

(V'|e-^^*|V') = J2 + ^ (e-("+^)*/2^ , (2) 

r,<a 

where are the resonances of H {i.e., the complex poles of {E — H)~^ in the lower half- 
plane) whereas Vj :— — 21m Ej are their widths and a and £ are suitable positive numbers. 
Herbst also showed that inf{rj} > 0, in such a way that Eq. (2) ensures the exponential 
decay of A{t; [ip]) as t oo. 

The purpose of the present paper is to investigate the decay law in a very simple (albeit 
non-trivial) model, namely, a one-electron atom in which the Coulomb attractive potential 
is replaced by a point-like attractive well — an idealization of a very short-range attractive 
interaction — and put under the infiuence of a static and uniform electric or gravitational 



^Another measure of decay, used for states initially confined inside a region (i.e., ip{x,t = 0) 
vanishes outside Ai), is the nonescape probability, defined as P{t) := Jjf^\ip{x,t)\'^ dx. Using 
Schwarz's inequality one can easily show that P{t) > \A{t)\'^, so that Khalfin's argument also rules 
out the exponential decay of P{t) for t — >^ oo. 
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field (for related simple models, see Rcfs. [9,10]). The model will be studied in D — 1,2,3 
space dimensions and will be shown to be exactly solvable in the one- and three-dimensional 
cases, whereas in the two-dimensional case it is solvable up to a quadrature. In spite of 
its simplicity, this model unravels some remarkable features that can be actually evaluated 
in detail and become worthwhile to be used as a paradigm with respect to more realistic 
situations, without any substantial change in the basic physical contents. To this concern, 
it is known that, in the absence of the uniform field, this model exhibits a bound state — 
see for instance Refs. [11-13]. 

It turns out that, once a background uniform field has been switched on, an infinite 
number of resonances arise in this model. In particular, the state vector that corresponds 
to the bound state in the absence of the uniform field is turned into a bona fide quasi-stable 
state for a sufficiently weak external field. For instance, if the bound state energy is of the 
order of 1 eV, the lifetime of the corresponding quasi-stable state in the presence of the 
Earth's gravitational field is much longer than the present age of the Universe; even in the 
presence of a rather strong laboratory static electric field, its lifetime is long in comparison 
to the typical time scales of atomic and condensed matter physics. On the other hand, 
very strong external fields are expected to create non-perturbative deviations from the naive 
exponential decay law. This has been observed in previous numerical studies of the present 
model [14-16] and will be qualitatively explained in this work. 

The paper is organized as follows. In Section II we first analyze the one-dimensional case, 
where a direct one-to-one correspondence takes place between the strength of the attractive 
potential well and the bound state energy. All the main features of the model are explicitly 
exhibited and discussed. In Section 111 we generalize our investigation to the two- and three- 
dimensional cases. Here the renormalization procedure is mandatory, in order to remove the 
ultraviolet divergences of the Green's functions. In so doing, the bound state energy in the 
zero-field case achieves a deeper physical meaning — it specifies the self-adjoint extension 
of the quantum Hamiltonian operator — whereas the renormalized coupling parameters 
become running auxiliary quantities. Mutatis m,utandis, all the main physical properties of 
the one-dimensional case are essentially recovered. In Section IV we draw our conclusions, 
whilst we defer some technical although important details to the Appendices. 

II. THE ONE-DIMENSIONAL MODEL 

Let us consider the Hamiltonian^ 

H = -— - X5(x) -Fx , A > , F > , (3) 

describing a particle interacting with an attractive (5-potential and a uniform background 
field. In the absence of the field (i.e., when F = 0) there is a single bound state 
with energy Eb — — the corresponding wave function being given by iPb{x) — 



^We use atomic units such that h = 2m = 1. 
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exp(— A|a;|/2). Once the uniform field is turned on, this bound state becomes un- 
stable, in the sense that A{t; [i^b]) -^Oast^ oo. The precise way in which this occurs will 
be the subject of this Section. 



A. Retarded Green's function 

The retarded Green's function G'^{E;x,x') is the solution to the differential equation 
{E-H)G+{E;x,x') = S{x-x') , EeC, (4) 

that satisfies the boundary condition 

lim G^{E;x,x') =0 for Im(E) > ; (5) 

>oo 

it is defined for Im(£') < by analytic continuation. The solution to Eq. (4) is known 
[12,16,17], but we shall derive it here for the sake of completeness. 

To solve Eq. (4), let us first consider the case A = 0; it can then be rewritten as 



(6) 



where 



p:=FV3(^^+|y (7) 

The solution to Eq. (6) that satisfies the boundary condition (5) is given by 

G+(p,p')=aAi(-p)%'-p) + 6CiW(-p)%-p') , (8) 

where Ai(a;) and Ci*'^-'(a;) := Bi{x) + iAi{x) are Airy functions [18] and 6{x) is the Heaviside 
step function. The coefficients a and b are fixed by the matching conditions at p = p': 

G+(p' + 0,p')=GoV-0,p') , (9) 

6>,G+(P,P')lp=.'+o - 9,G+(P,P')lp=p'-o = F-y' . (10) 

Solving these equations one finally arrives at 

Go+(p,p') = -7rF-V3 Ai(-p_) ci(+)(-p^) , (n) 

where 2p± := p + p' ± |p — p'|. 

In order to obtain G'^{E; x, x') for A 7^ 0, we rewrite Eq. (4) as an integral equation: 

/+00 
dy G^{E;x,y)XS{y) G^E;y,x') 
-00 

= G^{E] X, x') - A G^{E; x, 0) G+{E; 0, x') . (12) 
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Taking x = 0, solving for G~^{E; 0, x'), and reinserting the result into Eq. (12) yields the so 
called Krein's formula [11]: 

G*(E;xy)=GUE:.,x') - GUE^ rAS)G^(E:0...:') ^^^^ 



where 



y(A,£;)-i + G+(£;;0,0) . (14) 



B. Resonant-mode expansion of the propagator 

Prom G^{E; x, x') one can obtain the retarded propagator K'^{t; x, x') by a Fourier trans- 
formation: 

^ e-''''G+{E;x,x') . (15) 
-oo 27r 

It turns out that the following bound on G~^{E; x, x') holds true in the lower half-plane^ for 
\E\ sufficiently large (see Appendix A): 

\G^{E;x,x')\<C\E\-'/^e^^{\E\^'^{\x\ + \x'\)], \E\ ^ oo , (16) 

where C is a suitable constant. 

This bound allows one to close the contour of integration of (15) when t > with a semi- 
circle of infinite radius in the lower half-plane without changing the value of the integral. 
Using Cauchy's theorem, one then obtains the so-called resonant-mode expansion of the 
propagator [19,20]: 

K+{t- X, x') = Yl '^"(^) "fi-ix') , (17) 

n 

where the sum runs over the poles^ of G~^{E;x,x') located in the lower half-plane and the 
functions (pn{x) are given by 



[-a^G+(^;o,o)]V2 



(18) 

E—En 



^More precisely, the bound is valid only outside the sectors | arg(£') + 27r/3| < 5 and —5 < 
arg{E) < 0, with 6 > depending on As shown in Sect. II C and Appendix D, these regions 
contain poles of G'^{E;x,x') with arbitrarily large absolute values, where the inequality (16) is 
obviously false. One can, however, make 6 arbitrarily small by taking \E\ sufficiently large. 

^In writing (17) and (18) we have made use of the fact that the poles of G~^(E;x,x') are simple, 
as it can be explicitly checked by direct inspection. 
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The functions ipn{x) can be recognized as the so-called Gamow states [19,21]. On the one 
hand, just like the bona fide energy eigenfunctions, they satisfy the differential equation 
Hipn{x) = En(pn{x). On the other hand, the complex quantities do not correspond to 
the eigenvalues of the self-adjoint Hamiltonian operator and, moreover, the Gamow states 
are neither normalizable (not even in the sense of generalized functions, because they diverge 
when X oo) nor mutually orthogonal. 
Using Eq. (17) and the fact that 

/ + 00 
dx K^{t-x,x)ip{x\Q) , t>0, (19) 

-oo 

we can recast the non-decay amplitude A{t] := ('0|e"*^*|'^) in the form of a resonant- 
mode expansion: 

^(*;[V']) = ^^nC„e-^"*, (20) 

n 

where 

/ + 00 ^ /" + 00 

dxilj{x,Q) cpn{x), Cn-= / dxil)* {x,Ql) (pn{x). (21) 

-oo J —oo 

Notice that C„ ^ C* and \^n{x)\^ ~ exp(F-i/2r„ x^/^) as x ^ oo. It follows therefrom 
that the wave function ■?/^(a;, 0) of the initial state must decrease sufficiently fast at infinity in 
order that the coefficients C„, C„ exist. This condition is fulfilled by ■^(x, 0) = iPb{x)- This 
still leaves open the question of whether the series (20) converges. Here we shall assume 
that it does, at least in the Z^-topology. 



C. Poles of the Green's function 



The unperturbed Green's function Gq{E;x,x') is an holomorphic function of E, so that 
the poles of G^{E; x, x') are all given by the zeros of g{\ E). Inserting the explicit form of 
Gq{E; 0, 0) into Eq. (14) and noting that when F — there is a bound state with energy 
Eb — — we arrive at the following equation: 



Ai(-£)GiW(-£) = i-(-£^)-V2, 



:= EnF-^/^ . 



(22) 



For a given value of cb, Eq. (22) has an infinite number of solutions, all located in the lower 
half-plane. Some of them are shown in Fig. 1. They can be numbered according to their 
values in the limit eb — oo, which corresponds to a very weak field {F — > 0) or a very 
strong attractive interaction (A oo). One of the poles approaches the negative real axis 
and behaves asymptotically as (see Appendix B) 



£q ^ £b + iexp 



Sb -oo 



(23) 
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Its real part corresponds to the energy of the (unique) bound state of the atom in the 
absence of the uniform field. Its imaginary part is half the decay rate of the atom via 
tunneling through the potential barrier created by the external field. 

The other poles approach the zeros^ of Ai(— e), which are real and located on the positive 
real axis, 

lim En = -an , n e N , (24) 

£S-+-0O 

and of Ci^+\-e), 

lim = -an e"'^''^''^ , neN . (25) 

In Eq. (24), denotes the n-th zero of Ai{z); Eq. (25) follows from the identity Ci^~^\z) — 
2eW6Ai(;2e2^'^/3) [18]. 

If the external field is very weak, but nonvanishing — i.e., \eB\ ^ 1 — then the poles 
En with n > exhibit a small negative imaginary part (see Fig. 1). However, while Im(£o) 
approaches zero exponentially fast as Eb — —oo, one has Im(£„) ~ {—eb)~^ in the same 
limit (provided n is not very large, see Appendix C). This means that the transient effects 
associated to the poles with n > — and a fortiori those ones associated to with 
n < — disappear much faster than the corresponding effects associated to the resonance 

So- 

Looking at Fig. 1, one can notice that the imaginary parts of the first few poles En with 
n > have the same order of magnitude. This explains the short time oscillatory behaviour 
of \A{t; [V^])!^ observed in numerical studies [14-16] of the model (3) in the weak field regime: 
it is a consequence of the interference among the resonances associated with those poles. As 
a matter of fact, these resonances have a simple physical interpretation: when the external 
field F is turned on, it may excite the particle to a state of positive energy. Once excited, 
the particle is pushed to the positive x-direction by the field — recall that we are assuming 
F > — but it is scattered by the potential V{x) = —Xd{x). Because the potential is 
strongly attractive as A is very large, the transmission probability is small, so that the 
particle can bounce back and forth many times in the region to the left of the origin before 
it finally "jumps over" the potential well. 

Let us now examine the strong field regime le^l <^ 1. In this shown in Fig. 1, 

the decay rates {Ljl j G Z} form a monotonic decreasing sequence, with limj^oo^j = 0, as 
shown in Appendix D. Thus, in contrast with the class of potentials considered by Herbst 
[8], there does not exist a slower decaying resonance, which would eventually dominate the 
decay process. As a consequence, the decay is not asymptotically exponential: Va > 0, 
limt^oo \A{t; [ip])]"^ = oo. By the way, strictly speaking this result actually holds true 
even in the weak field regime \eb\ ^ 1, because lim„^oor„ = regardless the value of Eb 
(see Appendix D). In the weak field case, however, one should have to wait an extremely 
long time until a deviation from the exponential decay \A{t; [^])P ~ exp(— Lot) became 



^Note that the r.h.s. of Eq. (22) vanishes in the limit £b — — oo. 
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appreciable. Besides, \A{t; would be so small by then that such a deviation would be 
practically unobservable. 

The crossover from weak to strong field regime occurs at Eb ~ —1- At this value, 
Im(£o) ~ Ini(£:i) (see Fig. 1) — an indication that the two mechanisms of decay discussed 
above become equally important. 

III. THE TWO- AND THREE-DIMENSIONAL CASES 

A. Retarded Green's function 

We can use the same strategy employed in Sec. 11 A to solve the /^-dimensional version 
of Eq. (4) , which reads 

[E + V^ + A5(^) (x) + Fx]G+ {E; x, x') = ^^^^ (x - x') , (26) 

where x = (xi, . . . , xd) ■— {x, r) and E e C Thus we can formally write G^{E] x, x') as 
in Eq. (13), in which Gq (£^;x, x') denotes the solution to Eq. (26) in the case A = 0. The 
latter can be written as 

Gt{E-, X, x') = j e'>^-(^-') g^{E, k; x, x') , (27) 

where QQ{E,\i.]x,x') satisfies 

(^-\^ + ^ + Fx^ g+{E, k; X, x') = S{x - x') . (28) 

This has precisely the form of Eq. (4) with A = and £'—>£' — k^, the solution to which is 
given by Eq. (11). Inserting it into Eq. (27) we finally obtain 

G+(£^;x,xO = -vrF-V3| ^^^e'''-(-')Ai(-p_)CiW(-p+) , (29) 
where now p := F^^^ [x + {E — k^)/F] and 2p± p + p' ± \p — p'\. 

B. Renormalization 

In contrast with the one-dimensional case, the Green's function is ill-defined at coincident 
points for D > 2. Indeed, after setting x = x' = in Eq. (29) and performing the angular 
integration, we obtain 

POO 

G+{E; 0, 0) = -Cd / Ai{q) Ci^+\q) A;^-^ dk , (30) 

Jo 

where 
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D '■- 



22-D^(3-r')/2 

- l)/2] ' 



q := F-2/3(k2 - E) 



Since 



Ai(g)Ci(+)(g) 



-1/2 ^1/3 



27r 



for k 



oo 



(31) 



(32) 



the integral in Eq. (30) turns out to be ultraviolet divergent in D > 2. In the two- and 
three-dimensional cases the divergence can be absorbed through a redefinition of the coupling 
parameter A. To do this we follow the same procedure employed in [22]. Let us introduce 
a cutoff A in the upper limit of integration in Eq. (30) and add to the resulting expression 
the following integral: 



c 



D 



A V.D-2 



k^'-'dk 



27r Jo y^T^ ' 



(33) 



which contains the arbitrary momentum scale /x > 0. At the same time, we subtract /l)(A, h) 
from A^^ and define the rcnormalized coupling parameter as 



M/x)]-^:= hm [X-'-IniA,i,)], 

A— >oo 



(34) 



where it is understood that A depends on A in such a way that the limit exists. In this way, 
the denominator of the Krein's formula (13) is replaced by an expression that is finite when 
the cutoff is removed: 



lim goiX, E) = ^ - 



:= gD{XR,ii,E) . 



Ai(g) CiW(g) 



k^'-^dk 



(35) 



In the next two subsections we shall analyze this expression separately in D = 2 and 
D = 3 dimensions. It turns out that the latter is simpler than the former, so we discuss it 
first. 



C. Three-dimensional case 

In D = 3 the integral in Eq. (35) can be computed in closed form^ yielding (e :— EF~'^^^) 
1 n 1 



^3(A^, ^,E)^^-^-- F^l^ e Ai(-£) CiW (-£) + Ai'(-£) CiW'(-^) 



(36) 



Using the asymptotic expressions of the Airy functions for large argument [18], one can 
easily show that in the limit F — > the expression above is reduced to 



'/ 2/12/2 dx = xyiy2 — y'iy'2 for any two solutions of Airy's equation y" — xy = 0. 
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Thus, provided Aj? > ^'Kj ii-, the quantity QsIXr, has a real zero given by 

2 



(37) 



-Er — 



47r 



Aj?(//) 



(38) 



which can be identified as the energy of the unique bound state of the system. It is worthwhile 
to remark that the bound state energy is a physical quantity and turns out to be independent 
of the arbitrary scale /j,. From this physical requirement one can readily obtain the flow 
equation for the renormalized running couphng parameter: 



1 + (/x-/xo)[AR(/xo)/47r] 



(39) 



which exhibits asymptotic freedom, i.e., Afl(/^) — > as — > oo. 

After setting Eb '■— F~'^I^Eb we can rewrite the resonance equation g^{Xji, ^,E) = in 
the form 



-(-£B)'/' + £Ai(-£) CiW(-£)+Ai'(-£) Ci(+)'(-£) = , 

TT 



(40) 



that generalizes Eq. (22) to the three-dimensional case. As in the one-dimensional case, Eq. 
(40) has a solution Eq that tends asymptotically to £b in the weak-field regime (see Appendix 
B): 



£o ~ £b i 1 + -;{-^b) 



-3/2 



exp 



SB 



-OO . 



(41) 



It has the same physical interpretation of its one-dimensional counterpart — see Eq. (23). 

In addition to Eq, Eq. (40) has an infinite number of solutions. Some of them are shown 
in Fig. 2 for three different values oi Eb- Their distribution in the complex e-plane bears 
some resemblance with the one-dimensional case (see Fig. 1); in particular, they approach 
asymptotically the half-lines arg(£) = — 27r/3 and arg(£) = (see Appendix D). There are, 
however, two important differences: 

(i) for fixed n > we have that^ limeg_^„oo Im(en) ^ 0, as shown in Fig. 2, which clearly 
exhibits that the larger \eb\ the farther is e„ from the real axis; 

(ii) there is no clear distinction between the weak and strong field regimes — Fq is always 
smaller than Fi, even ioi Eb ^ 0. 

The first difference has a simple geometric interpretation: in three dimensions a particle 
can avoid a localized obstacle by going around it. Hence, the field F can easily detach a 
particle with positive energy from a localized potential. 



Note that this fact is not in conflict with Eq. (D6), which is vaUd under the condition that 
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D. Two-dimensional case 



Let us now finally discuss the two-dimensional case. To this concern, it is important to 
realize that the integral in Eq. (35) is no longer expressible in closed form when D — 2. Here 
we shall content ourselves with deriving an asymptotic expression for Eq in the limit F ^ 0. 
To this aim let us assume that |£o| is large and close to the negative real half-axis. In this 
case — see Eq. (31) — |g| is large and | arg(g)| ~ tt for all k G [0, oo), hence the following 
approximation is uniformly valid in the range of integration in Eq. (35) [18]: 



Ai{q) Ci(+)(g) 



fV3 r i 

, < 1 H — exp 



~F-' {k'-E) 



3/2 



(42) 



Inserting this into Eq. (35) we obtain, in = 2, 

1 1 / „2 



where 



HE) :-- 



A 



R 



dk 



^ln(-^ 
Att \ E 



Vk^-E 



exp 



{k' - E) 



3/2 



(43) 



(44) 



Consistently with our assumptions on E and F we can compute HE) using the saddle-point 
approximation and obtain 



HE) 



dk 



exp 



TlF 



£;)-3/4 



-If- 

3 



.£;)3/2 + ^ (_£;)i/2A;2 



exp 



F-^{-E) 



3/2 



(45) 



In the limit F — > 0, the integral I^iE) vanishes and g2{XR, fi, E) has a single real and 
negative zero Eb, corresponding to the energy of the bound state in the absence of the 
external field: 



Eb — — A*^ exp 



Att 



^'Ril^) 



(46) 



We notice that in D = 2 a bound state exists — provided, of course, F — — even if the 
renormalized strength of the point-like potential is negative, in which case one could naively 
expect the potential to be repulsive. 

As the bound state energy must be independent of the arbitrary scale ii, one can readily 
obtain the flow equation for the renormalized running coupling parameter, that now reads 



l + [Ai?(/^o)/27r]ln(/.//.o) ' 
leading again to asymptotic freedom. 



(47) 
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Now, let us consider Eq. (43) in the case of a weak field F. Using Eq. (46), we can 
rewrite it as 



g2{EB,E) — 



(48) 



An approximate solution to the equation (-Eb, -E) = is given by Eq — Eb [ 1 + ihiEs) ]; 
in terms of the dimensionless variable e — EF~'^/^ we obtain, cf. Eq. (45), 



-3/4 



exp 



£b 



-00 . 



(49) 



Moreover, as in the one- and three-dimensional cases, it is possible to show that an 
infinite number of resonances arise as solutions to the equation g2{^R- P - E) = 0, approaching 
asymptotically the half-lines a.Tg{E) = and a.Tg{E) = — 27r/3 when \E\ 00. 



IV. CONCLUSIONS 

In this paper we have analyzed the ionization of a very simple (though non-trivial) model 

atom submitted to the influence of a uniform static field. The model we have considered here 
is that of a one-electron atom in which the Coulomb interaction between the electron and 
the nucleus is replaced by an attractive short range (in fact, point-like) interaction. We have 
analyzed the problem in D = 1, 2 and 3 spatial dimensions. In spite of its simplicity and of 
the fact that — due to the external field — the Hamiltonian is not bounded from below, the 
study of the present model is far from being academic as it allows to grasp the basic features 
of the quantum dynamical behaviour of many realistic physical systems. In particular, its 
main prediction is a sensible deviation, in the strong field regime, from the naively expected 
exponential decay law of the survival probability of the bound state after the external field 
is turned on. Actually, more or less important deviations from the exponential decay law do 
occur even when the field is weak, specially in its short-time behaviour, with the presence 
of oscillatory transient effects (which are more pronounced in D = 1). Such deviations are 
caused by the presence of a purely continuous spectrum and the appearance of an infinite 
number of resonances once the uniform field is switched on. Deviations from the exponential 
decay law are also expected for very large times; this, however, may be an artifact of the 
model studied here, since for more realistic potentials one can prove asymptotic exponential 
decay [8]. (On the other hand, as noted before, the survival probability would be so small 
when such deviations took place that they would be practically unobservable.) 

An important development of the present investigation, which will be presented else- 
where, is the generalization of our analysis to the additional presence of a uniform magnetic 
field. In this way, it might be eventually possible to precisely evaluate the lifetimes of the 
so called non-conducting states — within the Integer Quantum Hall Effect (IQHE) conven- 
tional terminology — and to explicitly verify the widely popular picture according to which 
the presence of impurities, described in the simplest way by point-like attractive wells, gives 
rise to the plateaux formation in the IQHE [23] . 
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APPENDIX A 



In this Appendix we shall sketch the proof of the bound (16) on G~^{E]x,x) in D = 1. 
Let us first examine the asymptotic behavior of G~^{E; x, x') ioi \E\ —>■ oo in the sector — tt < 
arg(£') < — 27r/3. Since arg(p) — > arg(£') as oo, cfr. Eq. (7), we have |arg(— p)| = 

I arg(e*'^p)| < tt/S for large enough, so that [18] 



Ai(-p)~l7r-V2(_p)-i/4exp 



(Al) 



CiW(-p)-^7r-^^'(-p)"'^'iexp 



+ -exp 



ii-pr 



(A2) 



Inserting Eqs. (Al) and (A2) into Eq. (11), and dropping the second term in curly brackets 
in Eq. (A2) since it is negligible compared to the first one, we obtain 



(?+(£;; x,x')~^F-V3(p_p^)-i/4exp 
In addition, taking Eq. (7) into account, we have 



2i 



3/2 



-pT) 



(A3) 



p^/^r.F-'E^/^ + -E^/^x , 
2 



^ oo , 



so that 



(A4) 



G+{E;x,x') ~ ^£;-V2 exp {-iE'/^\x - x'\) . 



(A5) 



Inserting this expression into Eq. (13) and using the fact that lni{E^^'^) < 0, one can easily 
show that there is a positive constant C such that \G~^{E;x,x')\ < C\E\^^^'^ for \E\ suffi- 
ciently large and — tt < arg(£') < — 27r/3. Notice that Eq. (16) is a trivial consequence of 
this inequality. 

Let us now examine the asymptotic behaviour of G^{E;x,x') in the sector — 27r/3 < 
arg(£^) < 0. For this purpose it is convenient to rewrite Eq. (13) as 



G+{E;x,x') = 



G^{E;x,x') + XR{E;x,x') 
1 + \G+{E; 0, 0) 



(A6) 



where 
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R{E; X, x') := G^{E; 0, 0) G^{E; x, x') - G^{E; x, 0) G^{E; 0, x') . 



(A7) 



Again, since arg(p) — arg(£') as \E\ — oo, we have | arg(p)| < 2n/3 for \E\ large enough, 
in which case we have [18] 



Ai(-p)~7r-VV^/^sinQ//^ + ^) , 



(A8) 



exp 



SO that 



GUE;x,x')r.-F-'/'' (P-P+)-^/' 



X Ki exp 



2i 



3/2 , , 
P+ +P 



,3/2^ 



— exp 



2? / 3/2 3/2 



Using (A4) and neglecting the second term in curly brackets we obtain 



G^{E;x,x')r.-lE-'/'e^p 



-E^/^F-' + E'/\x + x') 
3 



(A9) 



(AlO) 



(All) 



The asymptotic behavior of R{E;x,x') depends upon the signs of x and x' . Let us 
first consider the case x,x' < 0. If we insert Eq. (11) into Eq. (A7), and use the identity 
f{p) f{p') = f{P-) we obtain {e := EF-'/^) 



R{E; X, x') = n^F-^/^ Ai(-p_) Ci^ (-e) 



X 



Ai(-£) CiW(-p+) - Ai(-p+) CiW(-£) 



(A12) 



Furthemore, inserting Eqs. (AS) and (A9) into Eq. (A12) and using Eq. (A4) we obtain 



R{E;x,x') — exp 



, , 1e^/^F-' + E'/''x. 



sm{E^/^x+) , x,x' <0 , (A13) 



where subdominant terms have been dropped. A similar analysis shows that 



R{E; X, x') ~ ^ exp 



2e'/^F-' + E^/^x+ 



sm{E^/^x^) , x,x' >0 , (A14) 



and that R{E;x,x') = for x < < a;' or < < x. Replacing Gq{E]x,x') and 
R{E; X, x') in Eq. (A6) with their asymptotic expressions, and using the inequalities | sin z| < 
exp(|Im^|), |Im£^^/^| < and \x ±x'\ < \x\ + \x'\, one can easily derive the bound 

(16). 
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APPENDIX B 



In this Appendix we derive the asymptotic expression of Eq in the weak field hmit. Let 
us first consider the one-dimensional case. If we assume that |£| » 1 and arg(£) a; — tt, then 
we may use the asymptotic expressions (Al) and (A2) for the Airy functions Ai and Ci^"*"^. 
Equation (22) then becomes 



-1/2 



1 + - exp 



-Sb) 



-1/2 



(Bl) 



This equation can be solved itcrativcly. As a first approximation, one may neglect the second 
term in square brackets, thus obtaining Eq^ Eb- In order to obtain the imaginary part of £o 
one must iterate once more: replacing e in the exponential with Eb and solving the resulting 
equation, one finds 



£o ~ £b i 1 +iexp 



£b 



-oo, D^l 



(B2) 



One can derive a systematic expansion in powers of Eb if one includes more and more terms 
in the asymptotic expansion of the Airy functions. In particular, the real part of the resulting 
expansion for Eq — F'^/^Eq agrees with the Rayleigh-Schrodinger perturbation series for the 
bound state energy when the external field is treated as a perturbation [14,16]. 

Following the same strategy, we can approximate Eq. (40) of the three-dimensional case 

by 



8£ 



exp 



. 



(B3) 



We can obtain an approximate solution to this equation using the iterative method employed 
above. This way, we finally arrive at the result displayed in Eq. (41). The result for the 
two-dimensional case is worked out in Subsection III D and is given in Eq. (49) . 



APPENDIX C 

In this Appendix we derive the asymptotic behavior of Ei in the weak field limit, Eb ^ 
— oo. As discussed in Subsection II C, in that limit the r.h.s. of Eq. (22) vanishes so that, 
to the lowest order, one has £i ^ — oi, where oi = —2.33810 ... is the smallest (in absolute 
value) zero of Ai(^). In order to obtain a more refined approximation, valid for a finite 
though large value of \eb\, we expand the l.h.s. of Eq. (22) in powers oi x = e + and, 
assuming that \x\ -C 1, we truncate the series and solve the resulting polynomial equation in 
X. The first non-trivial correction to the imaginary part of Ei is obtained when one truncates 
the series at 0{x^). In so doing, Eq. (22) is then approximated by a quadratic equation in 

ax^ + bx + c = Q , (CI) 
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where a = Ai'(ai) Ci(+)'(ai), ^ = -Ai'(ai)Bi(ai)_, and c = -(l/27r) (-£b)-^/2 < 1. Of the 
two solutions to Eq. (CI), x± = {—b d= Vfc^ -4ac)/2a, the one with the minus sign must be 
discarded, as it violates the condition that x ^ as Eb ^ — oo (c ^ 0). Expanding x+ in 
powers of c, we obtain 



Substituting a, 6 and c with their explicit expressions, we finally obtain 

1 (-^b)-^/^ 



(C2) 



Ai'(ai) Bi(ai) 27r 
Ci(+)'(ai) {-es)-^ 



Ai'(ai)2Bi(ai)3 47r2 



SB 



-OO 



(C3) 



An important consequence of this result is that lm(£:i) ~ [—SbY^ for e:^ — > — oo. 

One could be tempted to apply the reasoning above to any £„, n e N. However, there is 
an important caveat: the r.h.s. of Eq. (C2) is a good approximation to x+ only if |ac/6^| <C 1, 
or 



Ci(+)'(an) 



Ai'(a„) Bi(a„)^ 



<C |c| = 27r(-£B) 



1/2 



(C4) 



Using the asymptotic expressions of the Airy functions and of an — the n-th zero of Ai(2;) 



[18] — one can show that |a/&^| ~ 7r( 



vl/2 



. Hence, Eq. (C3) is also valid for n > 1 



(with the obvious substitution a\ an), provided |a„| -C \sb\- (See Appendix D for the 
asymptotic behavior of when \an\ ^ 



APPENDIX D 

In this Appendix we derive the asymptotic behavior of the resonances n 7^ 0, which 
are located very far from the origin in the complex e-planc. Let us first discuss the one- 
dimcnsional case. Assuming that \e\ ^ 1 and 6 := arg(£) ^ 0, we are allowed to use the 
asymptotic expressions (A8) and (A9) to the aim of approximating Eq. (22) by 



-1/2 



exp 



■ IS 



3/2 



+ i 



-1/2 



(Dl) 



If we further assume that \e\ ^ |£s|, we may neglect the second term in square brackets; 
the resulting complex equation is then equivalent to the following pair of real equations: 



4 , ,3/2 . 3^ 1 
- \e\ ' sm — — - In 
3 ' ' 2 2 



SB 



4 , 
3'^ 



<%io 3^ 9 

Y'^ cos — - - + 2n7r , 



n e N . 



(D2) 



Assuming n ^ 1 and \e\ large, one can easily solve these equations, obtaining 
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4 ' ' 



e ~ Sr. 



(D3) 



Since 16*1 -C 1, we can write e = \6\ e'^ ~ \e\ (1 + i9), so that 

n > 1, D = 1 . 



£n ~ Sn - ^-S^^/^ In 



(D4) 



Next we consider the case |£| ^ 1 and 9 :— arg(£) — 27r/3. Using the very same approxi- 
mations we readily come to the following estimate 



-2i7r/3 



+ T s-'/' In 



n > 1, D = 1 . 



(D5) 



A straightforward generalization of the above treatments to the basic resonance equation 
(40) in the three-dimensional case eventually leads to the following asymptotic expressions: 



£n ~ - 2 «n'/' 1^ (44/^) , £_„ ~ e-^^/^ 8*^ ; 71 » 1, L> = 3 . (D6) 
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FIG. 1. Poles of the Green's function in the complex e-plane (e_3 to eg, clockwise) in the 
one-dimensional case: eb = —10 (o), eb = —1 (O), £b = —0.1 (□), and £b = —0.01 (+). The 
dashed line corresponds to the half-line arg(£) = — 27r/3. (Angles appear distorded in this plot 
because the real and imaginary axes have different scales.) 
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FIG. 2. Poles of the Green's function in the complex £-plane (eo to £10, from left to rig 
the three-dimensional case: eb = —10 (o), = —1 (o), and eb = —0.1 (□). 
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